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Abstract 

The maximum likelihood method offers a standard way to estimate 
the three parameters of a generalized extreme value (GEV) distribu- 
tion. Combined with the block maxima method, it is often used in 
practice to assess the extreme value index and normalization constants 
of a distribution satisfying a first order extreme value condition, as- 
suming implicitely that the block maxima are exactly GEV distributed. 
This is unsatisfactory since the GEV distribution is a good approxi- 
mation of the block maxima distribution only for blocks of large size. 
The purpose of this paper is to provide a theoretical basis for this 
methodology. Under a first order extreme value condition only, we 
prove the existence and consistency of the maximum likelihood esti- 
mators for the extreme value index and normalization constants within 
the framework of the block maxima method. 

Key words: extreme value index, maximum likelihood estimator, 
block maxima method, consistency. 
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1 Introduction and results 

Estimation of the extreme value index is a central problem in extreme 
value theory. A variety of estimators are available in the literature, 
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for example among others, the Hill estimator [8], the Pickand's esti- 
mator |13| . the probability weighted moment estimator introduced by 
Hosking et al. |10| or the moment estimator suggested by Dekkers et 
al. [5]. The monographs by Embrechts et al. [7], Beirlant et al. [2] 
or de Haan and Ferreira [1] provide good reviews on this estimation 
problem. 

In this paper, we are interested on estimators based on the maxi- 
mum likelihood method. Two different types of maximum likelihood 
estimators (MLEs) have been introduced, based on the peak over 
threshold method and block maxima method respectively. The peak 
over threshold method relies on the fact that, under the extreme value 
condition, exceedances over high threshold converge to a generalized 
Pareto distribution (GPD) (see Balkema and de Haan [1]). A MLE 
within the GPD model has been proposed by Smith [18J. Its theoreti- 
cal properties under the extreme value condition are quite difficult to 
analyze due to the absence of an explicit expression of the likelihood 
equations: existence and consistency have been proven by Zhou [21] . 
asymptotic normality by Drees et al. [6]. The block maxima method 
relies on the approximation of the maxima distribution by a general- 
ized extreme value (GEV) distribution. Computational issues for ML 
estimation within the GEV model have been considered by Prescott 
and Walden |14^ [To] , Hosking [9 J and Macleod [UJ . Since the support 
of the GEV distribution depends on the unknown extreme value index 
7, the usual regularity conditions ensuring good asymptotic properties 
are not satisfied. This problem is studied by Smith |17] : asymptotic 
normality is proven for 7 > —1/2 and consistency for 7 > — 1. 

It should be stressed that the block maxima method is based on the 
assumption that the observations come from a distribution satisfying 
the extreme value condition so that the maximum of a large numbers 
of observations follows approximatively a generalized extreme value 
(GEV) distribution. On the contrary, the properties of the maximum 
likelihood relies implicitely on the assumption that the block max- 
ima have exactly a GEV distribution. In many situations, this strong 
assumption is unsatisfactory and we shall only suppose that the un- 
derlying distribution is in the domain of attraction of an extreme value 
distribution. This is the purpose of the present paper to justify the 
maximum likelihood method for the block maxima method under an 
extreme value condition only. 

We first recall some basic notions of univariate extreme value the- 
ory. The extreme value distribution distribution with index 7 is noted 
G 7 and has distribution function 

F 7 (x) = exp(-(l + 7 j;)- 1/7 ), 1 + 72; > 0. 



2 



We say that a distribution function F satisfies the extreme value con- 
dition with index 7, or equivalently that F belongs to the domain of 
attraction of G 7 if there exist constants a m > and b m such that 

lim F m (a m x + b m ) = FJx), xeR. (1) 

That is commonly denoted F € D{G^). The necessary and sufficient 
conditions for F £ D{G~ f ) can be presented in different ways, see e.g. 
de Haan [3j or de Haan and Ferreira [H chapter 1]. We remind the 
following simple criterion and choice of normalization constants. 

Theorem 1. Let U = (^j-rp^j be the left continuous inverse function 
of 1/(1 — F). Then F € D(G~) if and only if there exists a function 
a(t) > such that 

U(tx)-U(t) xT-1 
lim — - — —— — — = , for all x > 0. 

t-s>+oo a(t) 7 

Then, a possible choice for the function a(t) is given by 

( 7 U(t), 7 > 0, 

a(t) = \ -j(U(oo)-U(t)), 7 <0, 
[ U(t) - t- 1 U(s)ds, 7 = 0, 

and a possible choice for the normalization constants in (|TJ) is 

a m = a(m) and b m = U{m). 

In the sequel, we will always use the normalization constants (a m ) and 
(b m ) given in Theorem [TJ Note that they are unique up to asymptotic 
equivalence in the following sense: if (a' m ) and (b' m ) are such that 
F m (a' m x + b' m )^ Ey(x) for all then 

lim < = i and i im b' m -b m = Q 

rn— 5>+oo a m m— i>+oo a m 

The log-likelihood of the extreme value distribution G 7 is given by 
£ 7 (x) = -(1 + 1/7) log(l + jx) - (1 + ~/x)- lh , 

if 1 + jx > and —00 otherwise. For 7 = 0, the formula is inter- 
preted as £q(x) = —x — exp(— x). The three parameter extreme value 
distribution with shape 7, location \i and scale a > has distribution 
function x 1— > F^{ax + fi). The corresponding log-likelihood is 



(X yfj, \ 

— — J - log a. 



The set-up of the block maxima method is the following. We con- 
sider independent and identically distributed (i.i.d.) random variables 
(Xj)j>i with common distribution function F E D(G l0 ) and corre- 
sponding normalization sequences (a m ) and (b m ) as in Theorem [TJ 
We divide the sequence (Xj)j>i into blocks of length m > 1 and define 
the fe-th block maximum by 

M k , m = max(X (fc _ 1)m+1 , . . . ,X km ), k > 1. 

For fixed m > 1, the variables (Mfc m )fc>i are i.i.d. with distribution 
function i 7 " 1 and 

=> G 7o asm-) +oo. (o) 

Equation ([3|) suggests that the distribution of M k ^ m is approximately a 
GEV distribution with parameters (70, b m , a m ) and this is standard to 
estimate these parameters by the maximum likelihood method. The 
logdikelihood of the n-sample (Mi m , . . . , M n m ) is 



1 ™ 

L n (7,a,/x) = -^£ (7)|t)(T )(M fc(m ). 



n 
fc=i 



In general, L n has no global maximum, leading us to the following weak 
notion: we say that (7n;M n ,cr n ) is a MLE if L n has a focaZ maximum 
at (7n;Mn ) 5 : n)- Clearly, a MLE solves the likelihood equations 

™" = ° ™ th (4) 

Conversely, any solution of the likelihood equations with a definite 
negative Hessian matrix is a MLE. 

For the purpose of asymptotic, we let the length of the blocks 
m = m(n) depend on the sample size n. Our main result is the 
following theorem, stating the existence of consistent MLEs. 

Theorem 2. Suppose F € D(G^ ) with 70 > —1 and assume that 

m(n) 

lim ~ = +00. (5) 
- logra 



-00 



Then there exists a sequence of estimators (jmV-m&n) an d a random 
integer N > 1 such that 

P[(7n, fin, oy,) is a MLE for all n > N] = 1 (6) 

and 

^ a.s. /J"n b m a.s. „ , 0~n a.s. .. ,_v 

7„ — ► 70, ► I) ana > 1 as n — > +00. (7j 
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The condition 70 > —1 is natural and agrees with Smith |17j : it 
is easy to see that the likelihood equation (j3J) has no solution with 
7 < 1 so that no consistent MLE exists when 70 < —1 (see Remark [3] 
below). Condition © states that the block length m(n) grows faster 
than logarithmically in the sample size n, which is not very restrictive. 
Let us mention a few further remarks on this condition. 

Remark 1. A control of the block size is needed, as the following 
simple example shows. Consider a distribution F £ D(G^ Q ) with 70 > 
and such that the left endpoint of F is equal to —00. Then for 
each m > 1, the distribution of the m-th block maxima M^^m has 
left endpoint equal to —00 and there exist a sequence m(n) (growing 
slowly to +00) such that 

hm mm = — 00 almost surely. (8) 

KKn a m 

The log-likelihood L n {^f, a, y) is finite if and only if 

. A. . M k , m -n\ 
mm 1 + 7 > 0, 



l<fc<n 

so that any MLE (■%,, /^m <?n) must satisfy 

mm 1 + 7 n '— > 0. 

i<fc<n V a n J 

Using this observations, one shows easily that Equation ([5]) is an ob- 
struction for the consistency (J7J) of the MLE. Of course, this phe- 
nomenon can not happen under condition ([5]). 

Remark 2. It shall be stressed that condition ([S]) appears only in 
the proof of Lemma 0] below. One can prove that under stronger as- 
sumptions on the distribution F € D(G 7o ), condition ([5]). This is for 
example the case if F is a Pareto distribution function: one checks eas- 
ily that the proof of Lemma U goes through under the weaker condition 
lim n _j. +00 m(n) = +00. Hence Theorem [2] holds under this weaker con- 
dition in the Pareto case. In order to avoid technical conditions that 
are hard to check in practice when F is unknown, we do not develop 
this direction any further. 

The structure of the paper is the following. We gather in Section 2 
some preliminaries on properties of the GEV log-likelihood and of the 
empirical distribution associated to normalized block maxima. Section 
3 is devoted to the proof of Theorem [2J which relies on an adaptation 
of Wald's method for proving the consistency of M-estimators. Some 
technical computations (proof of LemmaU]) involving regular variations 
theory are postponed to an Appendix. 
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2 Preliminaries 



2.1 Properties of the GEV log-likelihood 

We gather in the following proposition some basic properties of the 
GEV log-likelihood. We note x~ and x+ the left and right end point 
of the domain £~, i.e. 

(x~,x+) = {x € R; 1 + 7X>0}. 

Clearly, it is equal to (—00,-1/7), R and (— 1/7, +00) when 7 < 0, 
7 = and 7 > respectively. 

Proposition 1. The function £ 7 is smooth on its domain. 
1- If 7 < — 1? ^7 * s stricly increasing on its domain and 



lim t*«(x) = — 00 lim £ 7 (x) 



+00 if "f < — 1 
i/7 = -l 



5. If^>—l,£yis increasing on (a; , x*] and decreasing on [x*, x+ ), 
where 

(1 + 7)^-1 

lim £ 7 (x) = lim + ^ 7 (x) = —00 

anc? ^ 7 reaches its maximum £ 7 (x*) = (1 + 7)(log(l + 7) — 1) 
uniquely. 

Remark 3. According to Proposition [TJ the log-likelihood £ 7 has no 
local maximum in the case 7 < — 1. This entails that the log- likelihood 
Equation ([!]) has no local maximum in (— 00,— 1] x R x (0, +00) and 
that any MLE (7 n) A*n> 5n) satisfies 7 n > —1. Hence, no consistent 
MLE does exist if 70 < — 1. The limit case 70 = —1 is more difficult 
to analyze and is disregarded in this paper. 

2.2 Normalized block maxima 

In view of Equation Q, we define the normalized block maxima 

MMk m — b m 
k,m= ! , k>l, 
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and the corresponding likelihood 

1 n _ 
L n (7,cj,/x) = - ^^( 7)At , CT )(M fcim ( n) ). 

fc=i 

It should be stressed that the normalization sequences (fl^) and ibm) 
are unknown so that the normalized block maxima M^ m and the like- 
lihood L n cannot be computed from the data only. However, they will 
be useful in our theoretical analysis since they have good asymptotic 
properties. The following simple observation will be useful. 

Lemma 1. i^jniV-ni^n) * s a MLE if and only if L n has a local maxi- 
mum at (%, (£t n - b m )/a m ,a n /a m ). 

Proof. The GEV likelihood satisfies the scaling property 

^,H,<t((x ~ b)/a) = £^ af i+b,aa)(x) + log a 

so that 

L n (j, fx, a) = L n (7, — , — J - log a m . 

Hence the local maximizers of L n and L n are in direct correspondence 
and the lemma follows. □ 



2.3 Empirical distributions 

The likelihood function L n can be seen as a functional of the empirical 
distribution defined by 

1 " 

Jl — ^ lvl k.m 
k=l 

where S x denotes the Dirac mass at point 16R. For any measurable 
/ : R — > [—00, +00), we note P n [/] the integral with respect to P n , i.e. 

1 n ~ 

Pn[/] = -^/(M fe , m ). 
H k=l 

With these notations, it holds 

L n (j,fi,a) = P n [^ (7)Ai)ff )]. 
The empirical process is defined by 

1 n 

F B (t) = P„((-oo, *]) = -£ l {M k „<*}» * G R - 

k=l 
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In the case of an i.i.d. sequence, the Glivenko-Cantelli Theorem states 
that the empirical process converges almost surely uniformly to the 
sample distribution function. According to the general theory of em- 
pirical processes (see e.g. Shorack and Wellner |16] Theorem 1, pl06), 
this result can be extended to triangular arrays of i.i.d. random vari- 
ables. Equation ([3]) entails the following result. 

Lemma 2. Suppose F E Z)(G 7o ) and \im n ^ +00 m(n) = +oo. Then, 

sup |F n (i) — F 10 (t)\ — ^> as n — > +oo. 

t(zR 

This entails the almost surely weak convergence P n G 70 , whence 

P»[/]"G 7 „[/] aS n^+oo 

for all bounded and continuous function / : R — > R. The following 
lemma dealing with more general functions will be useful. 

Lemma 3. Suppose F G D(G 7o ) and \\m n ^ +00 m(n) = +oo. Then, 
for all upper semi- continuous function f : R — > [—00, +00) bounded 
from above, 

limsupP n [/] < G JO [f] a.s.. 

n— >+oo 

Proof of Lemma\Q Let M be an upper bound for /. The function 
/ = M — f is non-negative and lower semicontinuous. Clearly, 

¥ n [f]=M-¥ n [f] and G 70 [f] = M- G J0 [/] , 

whence it is enough to prove that 

liminfP n [/] >G 70 [/] a.s.. 

n— i>+oo 

To see this, we use the relation 

P n [/] = / f(¥t(u))du. 
Jo 

where is the left-continuous inverse function 

F^ = inf{x G R; ¥ n (x) >u}, u G (0, 1). 

Lemma [2] together with the continuity of the distribution function 
Fy entail that almost surely, F^~(u) — > F^(u) for all u G (0, 1) as 
n —¥ +00. Using the fact that / is lower semi-continuous, we obtain 

liminf /(F£>)) > /(*£(«)) u G (0, 1). 



S 



On the other hand, according to Fatou's lemma, 

liminf / f(¥^(u))du> [ lim inf / (F£~ {u))du. 
Combining the two inequalities, we obtain 

liminf f(WZ(u)) > I f(F^(u))du, 



whence 



liminf P n [/] >G 70 [/] a.s.. 

n— y+oo 



□ 



The next lemma plays a crucial role in our proof of Theorem [2j Its 
proof is quite technical and is postponed to an appendix. 

Lemma 4. Suppose F £ D(G 1Q ) with 70 > —1 and assume condition 
(J5j) is satisfied. Then, 

lim F n [£ 10 ] = G J0 [£ 10 ] a.s.. (9) 

n— >-+oo 

It shall be stressed that Lemma 0] is the only part in the proof of 
Theorem [2] where condition © is needed (see Remark [2]) . 



3 Proof of Theorem E 

We introduce the short notation = ( — 1, +00) x R x (0,+oo). A 
generic point of is noted = (7,^,0"). 

The restriction L n : — > [—00, +00) is continuous, so that for any 
compact K C 0, L n is bounded and reaches its maximum on K. We 
can thus define 0^ = (7^ ,Ji^ ,&n) such that 

6% = argmax L n (0). (10) 

9&K 

The following proposition is the key in the proof of Theorem [2j 

Proposition 2. Let 0q = (70,0,1) and K C be a compact neigh- 
borhood of 0o. Under the assumptions of Theorem^ 

lim n = 0q a.s.. 

n— >+oo 

The proof of Proposition [2]relies on an adaptation of Wald's method 
for proving the consistency of M-estimators (see Wald [20J or van der 
Vaart [19] Theorem 5.14). The standard theory of M-estimation is 
designed for i.i.d. samples, while we have to deal with the triangular 
array {(Mk t m)i<k<n, n > 1}. We first state two lemmas. 
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Lemma 5. For all 9 £ Q, G l0 [£g] < G^ [£g ] and the equality holds if 
and only if 9 = 9q. 

Proof of Lemma [JJ The quantity G 7o [£g — £g] is the Kullback-Leibler 
divergence of th GEV distributions with parameters 9q an d 6 and is 
known to be non-negative (see van der Vaart |19) section 5.5). It 
vanishes if and only if the two distributions agree. This occurs if and 
only if 9 = 6q because the GEV model is identifiable. □ 

Lemma 6. For B C 0, define 

£b(x) = sup £g(x), x € M. 

Let 9 G and B(9,e) be the open ball in with center 9 and radius 
e > 0. Then, 

nmG 7Q [£ B{ g t£) ) = G JO [£g). 

Proof of Lemma Proposition Q] implies 

£g{x) = £~ ( ({x — fj>)/a) — log a < m 7 — loga. 

One deduce that if B is contained in (1,7] x [a, +oo) x R for some 
7 > — 1 and a > , then there exists M(j,a) such that 

£g(x) < M(7, a) for all 9 € B, x € R. 

Hence there exists M > such that function M—£ B ^g^ is non-negative 
for e small enough. The continuity of 9 i— >■ ^e( 2; ) on implies 

lim e -)(x) = £g(x) for all x € M. 

Then, Fatou's Lemma entails 

G 70 [liminf (M - £ B (e,e))] < liminf G 70 [M - £ Bm ] , 

whence we obtain 

limsupG 70 [£ i3(e )] < G 70 [4]. 

e^O 

On the other hand, 9 G B(9,e) implies G 70 [£ B {e,e)] ^ G^oM- We 
deduce 

lim G 70 [£ B (0 te )] = G yo [£g]. 

□ 
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Proof of Proposition [H In view of Lemmas [5] and [6l for each 9 G K 
such that 8 ^ 9o, there exists eg > such that 

Fix 5 > 0. The set A = {0 € iT; \\9 — 9q\\ > 5} is compact and 
is covered by the open balls {B(9,eg),9 G A}. Let Bi = B(9i,£g ), 
1 < i < p, be a finite subcover. Using the relation L n (9) = P re [fg], we 
see that 

supL„(#) < maxP n fe]. 
6»eA i<«<p 

The function is upper semi-continuous and bounded from above, 
so that Lemma [3] entails 

limsupPn^J < G 10 [9i] a.s., 

n— !-+oo 

whence 

limsupsupL n (#) < max G 7o [#j] < G^^o] a.s.. (11) 

n->+oo 6eA 1<«<P 
According to Lemma[U P n [^ ] ^> G 70 [^ 70 ], so that 

liminf supL n (0) > G 7o [0 o ] a.s.. (12) 

Since 9% realizes the maximum of L n over K, Equations (|lip and (fT2|) 
together entail that 9^ G K\A for large n. Equivalently, — #o|| < ^ 
for large n. Since 5 is arbitrary, this proves the convergence 9^ ^> 9q 
as n — > +00. □ 

Proof of Theorem \^ Let X C O be a compact neighborhood of #0 
as in Proposition [2] and define 0^ by Equation (|10p . We prove that 
Theorem [2] holds true with the sequence of estimators 

(%,V>n,(7n) = (in \ a mfin +t,"ra^), n > 1. 

According to Lemma [TJ (7™, /2 n , a n ) is a MLE if and only if L n has a 
local maximum at 0„ = (7^, ju„ , cr^ ). Since 6^ = argmax^g^ L n (9), 
this is the case as soon as L n lies in the interior set int(i^) of K. 
Proposition [2] implies the almost surely convergence 9^ 9q which 
is equivalent to Equation ([7|). Furthermore, since 9q G int(iT), this 
implies 9^ G int(-fC) for large n so that (7 n , /2 n , a n ) is a MLE for large 
n. This proves Equation ([6]). □ 
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Appendix: Proof of Lemma 3 

We will use the following criterion. 

Lemma 7. Suppose F £ D{G l0 ) and lim„_> +00 m(n) = +00. We 
note Y m = £ 70 (a~ 1 (Mi jm — b m )). If there exists a sequence (a n ) n >% 
and p > 2 such that 

^nP(|y m | > On) < +00 and supE[|y m j p l{|y m |< Qn }] < +00, 
n>l 

then Equation ([9]) holds true. 

Proof of lemma\^ We note [x = G 70 [£ 7() ] and we define 



n 

Y k) m = ^ { a m( M k,m ~ &m)) and S n = ^ Y k>m . 



k=l 



With these notations, ([9]) is equivalent to n 1 S n ^4' fi. We introduce 
the truncated variables 



n 

Y k>m = Yk,mX{\Y kim \<a n } an d &n = Y k , m . 



k=l 

Clearly, 

nS n ^S n ] < ¥[Y k>m ^Y k>m for some ke {l,...,n}] 
< nP[|y m | > a n ], 

so that Sn>i nF [l^ml > a n] < +00 entails ^jPfS'n 7^ 5 n ] < +00. 
By the Borel-Cantelli Lemma, this implies that the sequences (S n ) n >i 
and (S n ) n >i coincide eventually, whence n~ 1 S n -4' /i if and only if 
n _1 5 n -4-' /i. We now prove this last convergence. 

We first prove that E[Yi im ] —> [i. Indeed, by the continuous map- 
ping theorem, the weak convergence ([3|) implies Y\ m => £^ (Z) with 
Z ~ Gy . Since P[Yi m 7^ Yi m ] converges to as n — > +00, it also holds 
Y-i. m =>■ lry (Z). Together with the condition sup n >^ 
this entails E[Yi jTO ] -»• E[£^ (Z)} = p. 
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Next, Theorem 2.10 in Petrov |12| provides the upper bound 
E[\S n -E[S n ]\ p ] < C(p)n p / 2 E[\Y 1<rn -E[Y 1<m \ p ] 
for some constant C(p) > depending only on p. Equivalently, 

EOn- 1 ^ - /i„H < C(p)n-r/ 2 E[\Y hm - n n \ p \. 
with ji n = E[Y"i m }. Furthermore, 

E[|ii, m - Hn\ p ] < 2^ 1 (E[|y 1 , m n + |/x n |f) 

is uniformly bounded by some constant C > 0. By the Markov in- 
equality, for all e > 0, 

PfK 1 ^ - Mn| > e] < e-PEiln^Sn - ^ n \ p ] < £- p C{p)Cn- p l 2 . 

Since p > 2, it holds 

J^POn- 1 ^-^ >e] <+oo 

n>l 

and the Borel-Cantelli Lemma entails n~ 1 S n — fi n ^4' 0. Since fi n — > fi 7 
we deduce n _1 S n -4' /U which proves the Lemma. □ 

Proof of Lemma^ We prove that there exists a sequence (a n ) and 
p > 2 satisfying 

^na; p <+oo (13) 

n>l 

and 

supE[(|y m | Aa„f] < +oo. (14) 

n>l 

The Markov inequality yields 

n\y m \ >a n ] <a- p E[(\Y m \ A a n ) p ] 
so that Equations (|13p and (fT4"|) together entail 

^nF(|y m | > o n ) < +oo and supE[|y m | p l { |y m 

\<a n }\ < +00, 

n>l 

This shows that Equations (|13p and ()14p together imply the assump- 
tions of Lemma [7] and prove Lemma |4) 

We first evaluate the quantity £[(11^1 A a n ) p ] from Equation (|14D . 
Recall that Y m = ^ 70 ((Mi jm — b m )/a m ). It is well known that the ran- 
dom variable Xi with distribution function F has the same distribution 
as the random variable F*~(V), with V a uniform random variable on 
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(0, 1). We deduce that the random variable M\^ m = V^Xj has the 
same distribution as F^(V m ), with V m a random variable with distri- 
bution mv m ^ 1 l(o,i) [v)dv (this is the distribution of the maximum of 
to i.i.d. uniform random variables on [0, 1]). Hence, 

E[(|y m | A a n f] = ['(KdF^iv) - b m )/a m )\ A ocnfmv^dv. 
Jo 

The relations U(x) = F*~(l — 1/x) and b rn = U{m) together with the 
change of variable v = 1 — 1/ (mi) yield 

/OO -I 
/ (K(U m (x))\ Aa n ) p (l - — ) m - 1 x- 2 dx 

U{mx) — U (to) 



'l/m 

where 

U m (x) 



We now provide an upper bound for the integral and we use the fol- 
lowing estimates. There exists a constant c > such that 



K(y)\ < 



c (i + 7oyr 1/7 °, y<o, 
(1 + I/70) log(l + 7 y) + 1, y>0 



Note that U m (x) is positive for x > 1 and negative for x < 1. Further- 
more, for all x > l/m and to > 2, 

(1 - —) m ~ 1 < exp(-(m - l)/(mx)) < exp(-l/(2x)). 
nix 

Using these estimates, we obtain the following upper bound: for to > 
too (too to be precised later), 

E[(|y m | Aa„f] <h+h+h (15) 

with 

rm /m 
'l/m 



rmo/m 

h = <exp(-l/(2x))x" 2 dx, 

A/m 

/ 2 = T C P(l + 7 o^m(x))~ p/7o exp(-l/(2 a ;))^ 2 dx, 

Jmo/m 

/oo ____ 
((1 + I/70) log (1 + 7 of/ m (x)) + 1) P exp(-l/(2x))x- 2 dx. 

The integral I\ can be computed explicitely and 

h < 4aPexp(-TO/(2TO )). (16) 
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To estimate I2 and ^3, we need upper and lower bounds for U m (x) and 
we have to distinguish between the three cases 70 > 0, 70 £ (—1,0) 
and 7 = 0. 

Case 70 > 0: According to Theorem [TJ the function U is regularly 
varying at infinity with index 70 > and 

~ , . _ U(mx) — U(m) U(mx) 

1 + l0 U m (x) = 1 + 7o — — = -777 — 

a m U{m) 

We use then Potter's bound (see e.g. Proposition B.1.9 in [3]): for all 
e > 0, there exists uiq > 1 such that for m > too and rax > tuq 

U ( mx 1 

(1 — e)x 70 min(x e , x~ e ) < ——. — — < (1 + e)x yo max(x e , x~ £ ). 

U (m) 

We fix e G (0, 70) and choose mo accordingly. Using the lower Potter's 
bound to estimate I2 and the upper Potter's bound to estimate -Z3, we 
get 

I 2 < f c p ((l -e)x 70+e )" p/70 exp(-l/(2x))x- 2 dx 

< c?{l-s)- pho f 1 x- 2 - p - p£ ^°eM-l/(2x))dx, 
Jo 

and 

/oo 
((1 + I/70) log ((1 + e)x^ +£ ) + if exp(-l/(2x))x- 2 dx. 

These integrals are finite and this implies that I2 and ^3 are uniformly 
bounded for m > rriQ. From Equations (I15p and (|16p . we obtain 

E[(|y m | A a n ) p ] < 4< exp(-m/(2m )) + C, 

for some constant C > 0. Finally, we set exp(— m/(2mo)) = 1, i.e. 
a n = exp(m/(p2mo))- Equation (11411 is clearly satisfied and 

na~ p = exp[log n — raj (2mo)] = exp[— (m/ (2mo log n) — 1) log n]. 

We check easily that the condition lim n ^ +00 ™^ = +00 implies Equa- 
tion (HHD- 



Case 70 < 0: It follows from Theorem Q] that the function t 1— > 
U(oo) — U(t) is regularly varying at infinity with index 70 < and 
that 

~ . . U(mx) — U(m) U(oo) — U(mx) 

1 + loU m (x) = 1 + 70^ — = j T( \ ff, ' ■ 

a m U(oo) — U(m) 
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Then, the Potter's bounds become: for all e > 0, there exists mo > 1 
such that for m > mo and mx > mo 

n ™ , «■ „ U(oo) — U(mx) , , . , 

l-ek 70 min x e ,x" £ < TT \ \ -7--^ < (1 + e x 70 max x e , aT £ . 

(7 (00) — (7 (m) 

Using this, the proof is completed in the same way as in the case 70 > 
with straightforward modifications. 

Case 70 = 0: In this case, Theorem B.2.18 in |4j implies that for all 
e > 0, there exists mo > 1 such that for m > mo and mx > mo, 



U (mx) — U(m) , 

log x 



< emax(x e ,x £ ). 

Equivalently, for m > mo and mx > mo, 

log a; — £max(x e ,x~ e ) < U m (x) < logx + e max(x e , x~ £ ). 

Using the lower bound to estimate I2 and the upper bound to estimate 
I3, we obtain 

l _ 

<? exp(— pU m (x)) exp(— 1 / (2x))x~ 2 dx 



mo/m 



< d 



I exp(— p\ogx + pex 6 — l/(2x))x 2 dx, 
Jo 



and 



/oo ^ 
(U m (x) + l) p exp(-l/(2x))x- 2 (ix 

/oo 
(logx + ex £ + l) p exp(-l/(2x))x _2 dx. 

For e £ (0,1/p), the integrals appearing in the upper bounds are fi- 
nite and independent of m. This shows that I2 and I3 are uniformly 
bounded for m > mQ. The proof is then completed as in the case 
70 > 0. □ 
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